Bottleneck crossover between classical and quantum superfluid turbulence 
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We consider superfluid turbulence near absolute zero of temperature generated by classical means, 
e.g. towed grid or rotation but not by counterflow. We argue that such turbulence consists of a 
polarized tangle of mutually interacting vortex filaments with quantized vorticity. For this system 
we predict and describe a bottleneck accumulation of the energy spectrum at the classical-quantum 
crossover scale I. Demanding the same energy flux through scales, the value of the energy at the 
crossover scale should exceed the Kolmogorov-41 spectrum by a large factor ln 10 ^ 3 (£/ao) (£ is the 
mean intervortex distance and do is the vortex core radius) for the classical and quantum spectra to 
be matched in value. One of the important consequences of the bottleneck is that it causes the mean 
vortex line density to be considerably higher that based on K41 alone, and this should be taken into 
account in (re)interpretation of new (and old) experiments as well as in further theoretical studies. 



Introduction 



Turbulence in superfluid liquids, such as 4 He and 3 He 
at very low temperatures, is an intriguing physical prob- 
lem by itself because it comprises a system where the clas- 
sical physics gets gradually transformed into the quan- 
tum one during the energy cascade from large to small 
scales [1, 2]. Recently renewed broad interest to this 
subject has been motivated by an impressive progress in 
experimental techniques and new results, which has led 
in (at least) conceptual understanding of classical and 
quantum limits of the superfluid turbulence, see, e.g. [3- 
15]. Our paper in turn, attempts to shed light on physics 
of the superfluid turbulence behavior in the intermediate 
region near the classical-quantum crossover scale. We 
will see that transition of the turbulent energy cascade 
from the classical to the quantum scales is accompanied 
by transition from strong hydrodynamic to weak wave 
turbulence with a bottleneck stagnation at the crossover 
scale. 

Generally, the superfluid turbulence near zero tem- 
perature (for an introduction, see [1-3]) can be viewed 
as a tangle of quantized vortex lines. If turbulence is 
produced by classical means and not by a counterflow 
then at the scales much greater than the mean intervor- 
tex distance £, the vortex discreteness is unimportant, 
so that the superfluid turbulence has essentially a clas- 
sical character described by the Kolmogorov-41 (K41) 
approach [3]. As we will see below, vortex lines in K41 
state are polarized, i.e. tend to be co-directed and orga- 
nized in bundles. Since there is no viscosity or friction 
in a superfluid liquid near zero temperature, the classi- 
cal energy cascade proceeds down the spectrum to the 
scale of order £ without dissipation, where the vortex 
discreteness and quantization effects become important. 
Even though some negligible part of the energy is lost, 
for example, by radiation of phonons generated due to 
slow vortex motions and intermittent vortex reconnec- 
tions, the dominant part of the energy proceeds to cas- 
cade below the scale I by means of nonlinearly interacting 
Kelvin waves [3, 12, 13, 16], which were theoretically pre- 



dicted in the 19th century [17] and first experimentally 
observed by Hall [18]. We emphasize that the fact that 
the turbulence is produced by classical means is impor- 
tant here, because resulting polarization inhibits further 
vortex reconnections and prevents rapid fragmentation 
into vortex loops with sizes smaller than £. Thus, the 
main cascade carrier below scale £ will be Kelvin waves 
which are generated by both slow vortex filament motions 
and fast (but more rare and localized) vortex reconnec- 
tion events. Such reconnections produce sharp bends on 
the vortex lines and, therefore, generate a broader range 
of wavelengths than the slow vortex motions. However, 
the spectrum of the reconnection forcing decays with the 
wavenumber k sufficiently fast, and could effectively be 
thought as a large-scale source of Kelvin waves located 
at the crossover [15]. Traditionally, the K41 spectrum 
is assumed to maintain its shape all the way down to 
the crossover scale, which, due to such an assumption, is 
calculated based on the K41 spectrum [3]. 

In this paper, we demonstrate that, in contrast to the 
traditional viewpoint, the classical turbulent spectrum 
cannot be matched to its quantum counterpart at the 
same value of the energy flux because this flux requires 
much stronger levels of turbulence to be able to propa- 
gate through scales in sparse distributions of quantized 
vorticity. This leads to a bottleneck accumulation of the 
energy spectrum near the crossover scale which, in turn, 
significantly changes the position of the crossover, £, (see 
Fig. 1) and the relationship between the energy flux and 
vorticity, which have widely been used in interpretation 
of experimental results. Notice, that the phenomenon of 
bottleneck accumulation between two energy-flux spectra 
of different nature is not peculiar to the superfluid tur- 
bulence and may occur, for example, in the atmosphere, 
ocean and magnetosphere. 



I. 



POLARIZATION OF THE VORTEX TANGLE 



It is important that turbulence we consider in this pa- 
per is generated by classical means, e.g. by a towed grid 
[11] or by rotation [5-7], but not by a counterflow. In 
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FIG. 1: (Color online). The energy spectra Eh in the classi- 
cal, k < and quantum, k > l/£, ranges of scales. Two 
straight blue (dark gray) lines in the classical range indicate 

the pure K41 scaling E k oc k~ 5 ^ 3 , Eq. (12), and the pure 
thermodynamic scaling Eh oc k 2 . For the quantum range, the 
red (gray) solid line indicates the Kelvin wave cascade spec- 
trum, Eq. (11) (slope —7/5), whereas the green (light gray) 
dash-dotted line marks the spectrum corresponding to the 
non-cascading part of the vortex tangle energy (slope —1). 

the later case the vortex tangle would be unpolarized 
and neither we would expect K41 spectrum for the scales 
greater than £ (K41 is polarized, see below) nor would we 
expect Kelvin waves to be important for the scales below 
£ (reconnections would be more important, see [19]). On 
the other hand, polarization of the vortex tangle allows 
to shape large-scale vortex motions characteristic to the 
K41 cascade, and it also inhibits local reconnections and 
makes Kelvin waves a dominant vehicle for the turbulent 
energy cascade toward the small scales. Thus, let us con- 
sider the phenomenon of vortex polarization in greater 
detail. 

Intuitively, polarized vortex tangle can be viewed as 
a set of vortex bundles, so that in each bundle the vor- 
tex filaments have the same preferential direction. The 
simplest way to achieve such a polarization is to subject 
the system to an external rotation or shear. However, 
as we will see below, even isotropic and homogeneous 
turbulence can be, and often is, polarized. 

Let us formalize this picture by giving a mathematical 
definition of the vortex tangle polarization. Consider a 
circular disk or radius R with randomly selected position 
of its center and its orientation in the 3D space. The 
velocity circulation over the contour of this disc, T(i?), is 
obviously equal to the quantum circulation k (see Eq. (5)) 
multiplied by the difference between the number of vor- 
tices crossing the disk in the positive and the negative 
direction with respect to the normal to the disk, 

T(R) = k(N + - N_). (1) 

The totally unpolarized system is represented by a vortex 



tangle in which every vortex line consists of a chain of 
small uncorrelated segments (as in [19]). In this case 
the disc crossings would be completely random, and the 
mean value of T 2 would be determined from the central- 
limit theorem. Namely, if the sign of each crossing is 
completely random and statistically independent of all 
the other crossings then the total circulation T has zero 
mean and the standard deviation equal to the standard 
deviation for the circulation of an individual crossing k 2 
times the total number of terms in the sum (i.e. the 
number of crossings), 

(r 2 ) = K 2 (N+ + N-) ~ k 2 (R/£) 2 , (2a) 

where £ is the mean intervortex distance. We will say that 
this state has zero polarization, P = 0. Thus, the polar- 
ization P can be defined as a degree of deviation from 
this unpolarized state. For example, in the completely 
polarized system all vortex lines would be in perfectly 
aligned state, e.g. N- = and N + > 0, so that 

(r 2 ) = k 2 (nD ~ K 2 (i?/£) 4 , (2b) 

We will say that in this state P = 1. Let us now define 
polarization P by interpolating between these two limits. 
Namely, we will assume that the system is in a scaling 
state such that 

(r 2 > = k 2 {N 2 + ) ~ K 2 {Rjiy (3a) 

with some constant index a. Then for this state the 
polarization is defined as 

P = 0-/2-1. (3b) 

Note that in principle one can have a vortex system in 
which P < 0, e.g. an ordered grid structure composed of 
alternating positive and negative vortices. However, the 
alternating periodic structures are unstable and would 
quickly break up due to reconnections. 

Polarization of turbulent states with power law spec- 
tra is considered in Appendix B describing three differ- 
ent cases. For very steep spectra P = 1, for very shallow 
spectra (including the thermodynamic state) P = 0, and 
there is a window of intermediate spectra (including K41) 
for which P depends on the spectral slope and, there- 
fore, contains a nontrivial information about the turbu- 
lent scalings. For K41 turbulence we have 

(T 2 ) K41 ~e 2 / 3 R s / 3 . (4) 

In this case (r 2 ) can also be obtained from the di- 
mensional analysis. Thus, for K41 turbulence we have 
ct = 8/3 and polarization P = 1/3. 

Therefore, the vortex tangle associated with the K41 
cascade state is polarized. Note that in presence of bot- 
tleneck (described below) there will also be a contribu- 
tion of the thermalized part of the spectrum. However, 
this part is much less that the K41 contribution for large 
R/£. On the other hand, at scale R ~ £ (and obviously 
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FIG. 2: A sketch of typical reconnection of vortex lines in po- 
larized vortex tangles. A clash of two vortex bundles results 
in a localized reconnection region. Kelvin waves generated by 
the reconnections propagate away from the localized recon- 
nection region and spread in space along the vortex lines. 



for R < £) the notion of polarization becomes vague and 
useless, so one should not attempt to find P for these 
scales. 

Significant polarization associated with K41 cascade 
at large scales leads to grouping of the adjacent vortex 
lines into bundles with predominantly parallel orienta- 
tion which obviously inhibits reconnections and which 
selects Kelvin waves to be the dominant carrier of the 
downscale energy cascade. This picture is self-consistent 
because, as we will see later, only weak Kelvin waves are 
needed to carry the energy cascade of the same strength 
as in the large-scale K41 part. Associated with such weak 
waves small bending angles will not allow the adjacent 
(collinear) vortex lines to approach each other and recon- 
nect. On the other hand, it should be emphasized that 
the dominance of Kelvin waves over the vortex recon- 
nections still remains a hypothesis, even though a very 
plausible one. In this picture, the regions of reconnec- 
tions are intermittent and limited to locations where two 
vortex bundles clash, see Fig. 2. Kelvin waves generated 
by such localized reconnections will spread in space along 
the vortex lines into the vortex bundles. Therefore the 
resulting wave distributions will be much less localized 
in space than the reconnections. Note also that the re- 
connections in this picture do not lead to a creation of 
vortex loops of size I or less and, therefore, cannot trig- 
ger a cascade of further fragmentation of such loops, as it 
would be the case in unpolarized tangles in counterflow 
experiments [19]. Assuming that the large-scale dynam- 
ics of strongly polarized systems is similar to the classical 
flows described by Euler equations, one could imagine a 
classical prototype process in which reconnections will 
intermittently occur in locations of (yet to be proven to 
exist) singularities of the Euler equations. One can also 
see a clear analogy with reconnections of magnetic field 
lines in MHD with Alfven or whistler waves being similar 
to Kelvin waves. 



II. KINETICS OF INTERACTING KELVIN 
WAVES 

Let us describe the statistics of Kelvin waves on thin 
vortex filaments, and their role as a carrier of the en- 
ergy cascade at the scales less than the interline separa- 
tion I. Here we briefly overview the results of Kozik and 
Svistunov (2004) on this problem [12] (hereafter referred 
to as KS-04) with modifications and clarifications, par- 
ticularly keeping an explicit account for the logarithmic 
factors which will be important for the effects found in 
our work. The motion of the tangle of quantized vortex 
lines can be described by the Biot-Savart equation (BSE) 
[1, 2] for the evolving in time radius vector of the vor- 
tex line element s(£,i), depending on the arc lengths £ 
and time t. When the typical interline spacing I is large 
in the sense A = ln(£/ao) ^S> 1 (do is the vortex core 
radius), this equation can be simplified by the so-called 
local induction approximation (LIA) [20]. Both BSE and 
its LIA can be written in the Hamiltonian form [16]: 

inw = SH{w,w*}/dw* , 

where w(z, t) = x(z, t)+i y(z, t) with x and y being small 
distortions of the almost straight vortex line along the 
Cartesian z-axis. The BSE and LIA Hamiltonians are: 

bse k 2 f{l + Re[w'*(z 1 )w'{z 2 )]}dz 1 dz 2 . 
ti = — / , , (5a) 

47r J V( z i - z 2 y + \w{ Zl )-w{z 2 )\ 2 

H 1A = ^ jV{l + \W(z)\ 2 dz, (5b) 
k = 2-Kh/m , 

where primes denote the z-derivatives, k is the quantum 
of velocity circulation, and m is the particle mass. With- 
out the cut-off, the integral in H , Eq. (5a), would be 
logarithmically divergent with the dominant contribution 
given by the leading order expansion of the integrand in 
small Z\ — z 2 , that corresponds to H , Eq. (5b). 

It is well known that LIA represents a completely inte- 
grable system and it can be reduced to one-dimensional 
nonlinear Schrodinger (NLS) equation by Hasimo trans- 
formation [21]. However, it is the complete integrability 
of LIA that makes it insufficient for describing the en- 
ergy cascade and which makes it necessary to consider 
the next order corrections within the BSE model. 

Assuming that the Kelvin wave amplitudes are small 
with respect to their wavelengths, i.e. w' <C 1 (the 
self-consistency of this assumption is checked by an es- 
timate of the nonlinearity parameter, see Eq. (16)), 
we can expand the Hamiltonians (5) in powers of w' 2 : 
H = Hq + H 2 + Hi + Hq + . . . Next step is to consider 
a periodic system with the period length C (C — > oo 
to be taken later) and to use the Fourier representation 
w(z, t) = kT 1 / 2 ^2 k a(k, t) exp(ifcz) in terms of which the 
Hamiltonian equation takes the canonical form 

ida{k,t)/dt = dH{a,a*}/da(k,t)* 
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with a new Hamiltonian H{a, a*} — H{w, w*}/£ = 7i 2 + 
H4+H.6 + ... With dj = a(kj,t): 



H 2 = ^u k a(k)a*(k) , U A = ^ Ti 2 , 34 aia 2 a* 3 al , 
k 12,34 

^ e = 1q ^2 ^123,456 aia 2 a 3 ala* 5 al . (6) 

123,456 

Here Uk is the Kelvin wave frequency and interaction 
amplitudes Ti2,34, 1^123,456 are functions of k\ . . . ki and 
fci . . . ke, correspondingly. Summations over k\ . . . fc 4 in 
H4 and over k\ . . . ke in He are constrained by k\ + k 2 = 
fc 3 + fc 4 and fc 4 + k 2 + k 3 = fc 4 + k 5 + k 6 , correspondingly. 
One gets for functions in the Hamiltonians: 



BSE 
BSE 

1,2;3,4 

w BSE 

^1,2,3;4,5,6 



nA(k) fc 2 /47r , A(k) = ln(l/fca ) , (7a) 

fclfc 2 /c 3 /c 4 [A(£; of ) + ^l,2;3,4]/47T , 
9fcifc 2 fc3&4&5fc6[ A ( fc cf) + ^1,2,3;4,5,6]/327TK, 



LIA 2 LIA 

uj k = nAk /47T, T x 2 . 3 4 = kA kik 2 k 3 ki/A'K 1 

LIA 

^i,2,3;4,5,6 = 9 Akik 2 k 3 kik^ke/i2TTK ■ (7b) 

Here fc e f is a mean value of wave vectors in the game and 
all functions F are of the order of unity, they depend of 
the ratios of involved kj to k c {. 

It is well known that four-wave dynamics in one- 

BSE LIA 

dimensional case with dispersion laws Lo k or Lo k is ab- 
sent because the conservation laws of energy and momen- 
tum allow only trivial processes with k\ — fc 3 , k 2 = fc 4 , or 
k\ = fc 4 , k 2 — fc 3 . However, nontrivial six- wave scattering 
processes of 3 — > 3 type are allowed. For weakly inter- 
acting waves this dynamics can be described in terms of 
correlation functions (\a(k,t)\ 2 ) = LT 1 n(k,t), with the 
help of a classical six- wave kinetic equation [12], shown 
below for the continuous limit C — > 00 and rij = n(kj,t): 



jf\Wk,l,2;3,4,5\ 2 Ws, 4,5;fc, 1,2 — .A/fc,l,2;3,4,5] 



dn k 7r 
~df _ 12 J 

x5(uJ k +LUl+UJ 2 -UJ3-UJ A -U)c,) (8) 

x5(k+ki + k 2 — fc 3 — fc 4 — k 5 ) dk\dk 2 . . . dk§ , 

■A/i,2,3;4,5,6 = ™1™2»1 3 («4«5 + «4"6 + n 5 n 6 ) • 

Here H 7 is the /uiZ interaction amplitude, that includes 
the bare 6-wave amplitude W and 72 contributions of 
the 2nd order in the 4-wave amplitudes of the order of 
T kl 23 /to k . Notably, that LIA has infinitely many inte- 
grals of motion due to its complete integrability. These 
integrals totally preserve system from dynamical evolu- 
tion: dn(k,t)/dt = for any n(k) distribution. With 6- 
wave kinetic Eq. (8), this is possible only if W k 1 2 . 3 4 5 = 
on the resonant manifold k + k\ + k 2 = k 3 + k± + fc 5 , 
ujk +wi + uj 2 = u>3 +W4 + W5. This means that the leading 
contribution to W k 12-3451 proportional to A (that coin- 
cides with W k ^ 2 . 3 4 5 ) also vanishes due to cancelations 



of the leading contribution to W k 1 2 . 3 4 5 with that origi- 
nating from the perturbative terms. Remaining terms in 
W k 4 2 . 3 4 5 can be presented as follows 

BSE , 

Wfe,l,2;3,4,5 = kik 2 k 3 kihk e $l, 2 ,3;4,5,6/47rK, (9) 

where some dimensionless function <&i,2,3;4,5,6 is of the 
order of unity and depends only on mutual ratios of k- 
vectors k\ . . . k 6 . This estimate differs from Eq. (7a) for 

W k 1 2 . 3 4 5 by absence of the large factor A. 

The kinetic equation (8) written for a single vortex 
filament has a stationary solution [12] with a constant 
energy density (per unit length) flux e. We reformulate 
this "KS-04" spectrum to the 3D vortex tangle system 
in terms of the rate of energy density (per unit mass) in 
the 3D space, e — ejpi 2 (p is the fluid density): 

n k ~ (£ 2 e) 1/5 n 2/5 \k\- 17/5 , KS-04 spectrum. (10) 

It should be mentioned, though, that the theory of the 
cascade energy spectrum (KS-04) was derived with an 
assumption that vortex lines in the tangle are not recti- 
linear and non-interacting. In the present work, having in 
mind that reconnections are dominated by the mean in- 
tervortex distance, we silently assumed that the interac- 
tions and non-rectilinearity of vortex lines become unim- 
portant at small scales. 



III. WARM CASCADES IN HYDRODYNAMIC 
TURBULENCE 



The energy density per unit mass for Kelvin waves of 
small amplitude is 

£ = lJ uj k n k dk/2n = J £ k dk/2-K, 

where L ~ £~ 2 is the vortex line density per unit volume 
and £k is one-dimensional energy density in the fc-space. 
Together with Eq. (10) this gives 



£ k ^A( K 7 e/f) L/t> \k\ 



8\VB it! -7/5 



(11) 



Note that the parameters e and i in Eq. (11) are mutually 
dependent. Their relation follows from the expression for 
the mean vorticity in the system of quantum filaments, 
~ kL ~ k£~ 2 , where (\uj\) is dominated by the 
classical-quantum crossover scale and its estimate is usu- 
ally based on the K41 spectrum, 



£ 2 / 3 |fc|- 5 / 3 . 



(12) 



which gives 



(M) 2 ~ k 2 £ k "dk ~ £ 2 / 3 r 4 / 3 , or e^^li\ 

However, this estimate is rather unprecise because the 
K41 spectrum cannot be matched to the Kelvin wave 
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spectrum Eq. (11) at the crossover scale and, as it is 
explained below, there exists a bottleneck. But since the 
bottleneck is on the classical side of the spectral range, 
and the mean vorticity is still dominated by the crossover 
scale, one can find the correct relation between e and I 
based on Eq. (11) instead of K41. This gives: 

(|w|> 2 ~l/^fc|fc=V/~A(/t 7 e /^ 16 ) V5 , or e~K 3 /A 5 £ 4 . 

This estimate is different from the standard one based on 
K41 by a large factor of A 5 . 

Now, from Eqs. (11) and (12), one can find the ratio 
of quantum and classical (K41) spectra of turbulence at 
the crossover scale k ~ l/l: 



'-ye/ t-yt 



A 10 / 3 » 1 . 



(13) 



This ratio shows that quantum turbulence of Kelvin 
waves requires much higher level of energy (by factor 
A 10 / 3 ) in order to provide the same rate of the energy 
flux (and the same rate of the energy dissipation) than 
in the hydrodynamic turbulence of classical fluid. The 
main reason of that is the "rigidity" of the vortex fil- 
aments which is reflected by factor A k in Eq. (7a) in 
the Kelvin wave frequency. This contributes a factor of 
A 8 / 3 into the ratio (13). The remaining factor of A 2 / 3 
originates from the fact that any one-dimensional system 
of interacting Kelvin waves described by the Bio-Savart 
equation is close to the fully integrablc LIA system, in 
which the dynamics of wave amplitudes is absent. Thus, 
in order to have the same value of the energy flux and 
continuity of the spectrum at the crossover scale, there 
must be a bottleneck pile-up of the classical spectrum 
near this scale by the factor of A 10 / 3 , and this will be 
described by a "warm cascade" solution in what follows. 

As we explained, the energy flux carried by classical hy- 
drodynamic turbulence with K41 spectrum (12) cannot 
fully propagate through the crossover region. Therefore, 
hydrodynamic motions with larger scales (smaller wave- 
vectors) will increase their energy up to the level £\/i, re- 
quired for Kelvin waves to maintain the same energy flux. 
As the result, for k < l/£ the spectrum of hydrodynamic 
turbulence, £ k , will not have the K41 scale-invariant 

K41 

form E k given by Eq. (12). To get a qualitative un- 
derstanding of resulting bottleneck we will use so called 
"warm cascade" solutions found in [14]. These solutions 
follow from the Leith-67 differential model for the energy 
flux of hydrodynamical turbulence, 

HD 

where F k is the 3-dimensional spectrum of turbulence. 
Generic spectrum with a constant energy flux can be 
found as the solution to the equation e k = e: 



F k = 



24e 



ll|fc| n /2 \ np 



T \ 3 / 2 



2/3 



(15) 



The large k range describes a thermalized part of the 
spectrum with equipartition of energy characterized by 
an effective temperature T, namely, T/2 of energy per a 
degree of freedom, thus, F k — T/np and £ k = Tk 2 /np. 
At low k, Eq. (15) coincides with K41 spectrum, Eq. (12). 

This "warm cascade" solution describes reflection of 
K41 cascade and stagnation of the spectrum near the 
bottleneck scale which, in our case, corresponds to the 
classical-quantum crossover scale. To obtain the spec- 
trum in the classical range of scales, it remains now to 
find T by matching Eq. (15) with the value of the Kelvin 
wave spectrum at the crossover scale £ k ~ k 2 /£. This 
gives T/p ~ k 2 1 ~ (/s 11 /A B e) 1 / 4 . 

Obviously, the transition between the classical and the 
quantum regimes is not sharp and in reality we should 
expect that a gradual increase of the role of the self- 
induced wave-like motions of individual vortex lines with 
respect to the collective classical-eddy type of motions 
of the vortex bundles. Thus, the high-wavenumber part 
of the thermalized range is likely to be wave rather than 
eddy dominated. However, the energy spectrum for this 
part would still be of the same k 2 form corresponding 
to the thermal energy equipartition. This picture relies 
on the assumption (justified below) that the self-induced 
wave motions have small amplitudes and, therefore, do 
not lead to reconncctions. 

The resulting spectrum including both the classical, 
the quantum and the crossover parts, is shown on Fig. 1 
as a log- log plot. Important, that at k > 1/t in addition 
to the cascading energy associated with Kelvin waves, 
there is also energy associated with the tangle of vortex 
filaments (shown on Fig. 1 by a green dash-dotted line). 
The energy spectrum of this part ~ \k\~ l , which is sim- 
ply a spectrum associated with a singular distribution of 
vorticity along ID curves in the 3D space [10] and does 
not support a down-scale cascade of energy. The cascad- 
ing and non-cascading parts have similar energies at the 
crossover scale, that is the wave period and the ampli- 
tude are of the order of the characteristic time and size of 
evolving background filaments. In other words, the scales 
of the waves and of the vortex "carcass" are not separated 
enough to treat them as independent components. This 
justifies the matching of the classical spectrum at the 
crossover scale with the Kelvin wave part alone ignoring 
the "carcass" which is valid up to an order-one factor. 
This also justifies the way of connecting the "carcass" 
spacing I to the cascade rate e. 



IV. WEAKNESS OF TURBULENCE AT AND 
BELOW SCALE I 

In principle, turbulent fragmentation cascade into de- 
creasingly small vortex loops can be an alternative to 
Kelvin waves as mechanism of the energy transfer be- 
low scale I. Dimensionally, one can obtain spectrum 
cx |fc| _1 which corresponds to such a cascade [16, 22] 
(and which occidentally coincides with the non-cascading 
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"carcass" spectrum discussed above) . The probability of 
such small-scale reconnections depends on statistics of 
vortex orientations and can be estimated only in the sim- 
plest case of totally unpolarized vortex tangle by adopt- 
ing a model in which every vortex line consists of short- 
correlated pieces, see e.g. [19]. This model is relevant for 
turbulence produced by a thermal counterflow but not 
to the case of polarized turbulence produced by classical 
means. Definitely, in the case when, at the microscopic 
level, the vortex lines are preferably parallel (presumably 
this is the case for superfluid turbulence in the rotating 
tank [5-7]), the reconnection scenario [16, 22] is irrele- 
vant, as we assumed in our approach. 

This picture is supported by estimation of the nonlin- 
earity parameter through comparison of the nonlinear 
frequency shift Aujk with the frequency itself: 

= Aujk ?S.fc |jfc fc 1_ , , 

which is obtained using Eqs. (7a), (10) and our estimate 
e ~ k 3 /A 5 £ 4 . Now we see that at the cross-over scale the 
nonlinearity is small: 

£ 1/e ~l/A~l/ln(*/ao)«l . 

Correspondingly, characteristic values of the bending an- 
gle a associated with Kelvin waves are also small, 

a~ l/>/A< 1 . (17) 

Hence, the KS-04 weakly nonlinear spectrum should 
dominate the Svistunov-95 reconnection spectrum [16]. 
Indeed, the mean wave amplitude at scale I is <~ at, i.e. 
too small for the adjacent vortex lines to "touch" and re- 
connect. This proves self-consistency of the picture of the 
cascade carried by interacting Kelvin waves without re- 
connections, but we emphasize that this picture assumes 
polarization of the vortex system on which Kelvin waves 
propagate. 

Similarly, in the thermalizcd region of scales the 
mean bending angle can be estimated as a ~ \/£ ~ 
(k£) 3 /\/A <C 1. Thus, the self-induced vortex line mo- 
tions gradually arising from the eddy-like collective mo- 
tions in the thermalizcd part take form of weakly non- 
linear Kelvin waves. The nonlinearity of Kelvin waves 
grows with k in the thermalizcd part reaching its peak at 
the crossover scale and decreasing in the Kelvin cascade 
range. 

Summary and Discussion 

In this paper, we suggested the following route for the 
development of the energy cascade: K41 — > "Warm Up" 
— > KS-04 spectrum with a very pronounced Bottleneck 
effect. This scenario is relevant for polarized vortex sys- 
tems resulting from forcings of a classical type, e.g. a 
towed grid or rotation, but not relevant to unpolarized 



vortex tangles produced by thermal countcrflows. In our 
arguments we relied on the fact that polarization sup- 
pressed the reconnection-fragmentation cascade. Classi- 
cally produced K41 turbulence is indeed polarized. How- 
ever, its polarization is not perfect and at this time we 
cannot exclude that in some specific cases, the reconnec- 
tion dynamics can suppress the bottleneck accumulation 
of energy. 

In this paper we predicted that the bottleneck on the 
classical-quantum crossover scale amplifies the spectrum 
at this scale by a large factor of A 10 / 3 with respect to 
K41. Correspondingly, the corrected estimate for the 
crossover scale which takes this bottleneck into account 
is I ~ («; 3 /A 5 e) 1 / 4 , which is A 5 / 4 times smaller than the 
standard estimate based on K41. Yet another way to 
reformulate the same thing would be to say that the ef- 
fective viscosity v' is reduced by a large factor of A 5 , 
i.e. 

u'kk/A 5 (18) 

(see e.g. [3] for definition of v' and explanation of its 
meaning). 

A comment is due about locality of the transition be- 
tween the classical turbulence and the Kelvin wave cas- 
cade. Due to a sharp kink- like nature of the vortex recon- 
nections generating Kelvin waves (see Fig. 2), one might 
think that the energy is injected into the Kelvin wave cas- 
cade over a wide range of wavenumbers (associated with 
a Fourier analysis of the kink), and conclude that the en- 
ergy spectrum in the quantum region should differ from 
fc -7 / 5 . However, as it was shown in Ref. [23], the Fourier 
transform of the kink decays with k fast enough for the 
direct cascade scaling to dominate. In other words, the 
reconnection forcing appears to be more or less equiva- 
lent to a low-frequency forcing at the intervortex scale 

e. 

To describe the shape of the bottleneck spectrum, wc 
used the "warm cascade" solution previously obtained in 
[14] based on the Leith-67 differential model Eq. (14), 
as it is the simplest model which can provide a clear 
qualitative understanding of the bottleneck phenomenon. 
Clearly, the differential model Eq. (14) exaggerates lo- 
cality of the interactions of scales in real Navier-Stokcs 
turbulence where the main contribution to evolution of 
£ k originates from a wider range of comparable scales 
q ~ k. Some authors claim that extended interaction 
triads with q between k/A and Ak (with A <~ 10) are 
most important [24]. If so, the transient region between 
K41 and the thermodynamically equilibrium spectra can 
be wider than the one predicted by the differential ap- 
proximation Eq. (14). To account for this effect, one can 
use a more sophisticated turbulence closure based on an 
integral rather than the differential equation, e.g. one 
of the traditional closures such as the Direct Interaction 
Approximation (DIA) or "Eddy-Damped Quasi-Normal 
Markovian" (EDQNM), as it was done in Ref. [25], or 
even simpler closure, suggested in Appendix A. 
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In this paper, we did not consider the effects of the mu- 
tual friction between the normal and superfluid compo- 
nents, thereby restricting our consideration to low tem- 
peratures (e.g. below IK for 4 He). At higher temper- 
atures the dissipation of energy by the mutual friction 
can exceed the energy transfer to Kelvin waves, which 
would make our analysis and conclusions inapplicable. 
This seems to be the case, for example in experiments 
described in [11]. 

At lower temperatures there is a clear lack of exper- 
iments on turbulence generated by classical means. In 
this respect, one could mention the 3 He experiment on 
turbulence generated by a vibrating wire at Lancaster 
[4] the authors of which found the value of the effective 
viscosity v' = 0.2k which appears to be much greater 
than our prediction (18). On the other hand, to obtain 
this value the authors used an estimate for the integral 
(energy containing) scale to be equal to the thickness of 
the turbulent region, d = 1.5 mm which in our opinion 
may not be the case for the oscillating grid setup. Lack- 
ing direct measurements of d, we could get guidance from 
the oscillating grid experiments in classical fluids, see e.g. 
[28] where the following estimate for the integral scale is 



given, 



d = 0.25(5/M) 1 / 2 z, 



where M is the mesh size, S is the amplitude of oscilla- 
tions, z is the distance from the grid. Taking the Lan- 
caster parameters, M = 50 /im, S = 1 /im, z = 1.5 mm, 
we get d = 50 /im. Estimating v' with this value of the 
integral scale would give i/'~2x 10~ 4 k, which would be 
consistent with the small A -5 coefficient in (18). How- 
ever, it is not possible to be more conclusive one way or 
another without more direct measurements of the turbu- 
lent parameters in this case. 

A lucky exception appears to be a new 3 He experiment 
on rotation generated turbulence, in which the bottleneck 
phenomenon appears to be important in understanding 
the observed propagation speed of the turbulent-laminar 
interface, see [26] for detailed explanations. 

In the final stages of modifying our paper according 
to the referee comments, our attention was drawn to 
a new preprint [29] where an alternative picture of the 
crossover turbulence was presented with bottleneck be- 
ing prevented by reconnections. The authors argued that 
the self-induced part of the vortex line velocity becomes 
larger than the classical collectively produced velocity in 
the vortex bundle already at the scale r ~ A 1 ' 2 Zq S> lo- 
From this they concluded that the polarized vortex bun- 
dles move randomly with respect to each other which 
leads to their random reconnections. In this respect we 
would like to re-emphasize our view that the fast self- 
induced motions take form of rapidly-oscillating Kelvin 
waves rather than of a chaotic motion of vortex bundles. 
Moreover, as expressed in our estimate (17), these Kelvin 
waves must have very small bending angles (1/A 1 / 2 or 
less) in order for the six-wave Kelvin cascade to carry 
the same flux as in the K41 (large-scale) part of the spec- 



trum. Such small bending angles are insufficient for the 
neighboring lines within a particular bundle to approach 
each other and to reconnect neither in the thermalized 
nor in the cascade range of scales. Thus, reconnections 
are limited to rather small volumes in-between of collid- 
ing large-scale bundles. 

On the other hand, as we have already said in this 
paper, polarization of K41 turbulence is not perfect and 
we may expect reconnections to play some role which 
would lead to certain modifications of the bottleneck phe- 
nomenon described in this paper. Relative role of the re- 
connections versus the Kelvin wave cascades is also likely 
to depend on the particular way of turbulence excita- 
tion. In particular, we may expect further reinforcement 
of the polarization, and therefore stronger suppression of 
reconnections, in rotating systems and in systems with a 
strongly sheared mean flow. 
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APPENDIX A: SIMPLE TURBULENT CLOSURE 

Here we propose a new model which could be viewed 
as the simplest (minimal) integral closure to be used 
in the future for an improved description of the tran- 
sitional bottleneck region. The model comprises in writ- 
ing the collision term St{Ffe} in energy spectrum balance 
dF k /dt = St{F k } as follows: 



St{F fc } ~ J 



oo „2 



q 2 dqp 2 dpS(k + q + p) 
2nk 2 (k 2 + q 2 + p 2 ) 



(Al) 



x k [k F q F p + q F k F p + p F q F k ] /( 7fc + 7? + lp ) , 

where k, q, and p are one-dimensional vectors varying 
in the interval (— oo, +oo), and -f k = \/\k\ 5 F k represents 
eddy-turnover frequency. The model (Al) differs from 
EDQNM by replacement of d 3 qd 3 pS 3 (k + q + p) with 
3-dimensional vectors fc, q, and p by q 2 dqp 2 dp 5(k + 
q + p) I (k 2 + q 2 + p 2 ) with one-dimensional vectors, by 
a simpler form of j k , and by one-dimensional version of 
the interaction amplitude (V^ 7 => k). 

The model (Al) satisfies all general closure require- 
ments: it conserves energy, J k 2 St{F k }dk = for any 
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F k , and St{-Ffc} = on the thermodynamic equilib- 
rium spectrum F k =const and on the cascade spectrum 
Fk oc |/c| -11 / 3 . Importantly, the integrand in Eq. (Al) 
has the correct asymptotical behavior at the limits of 
small and large q/k as in the sweeping- free Bclinicher- 
L'vov representation, see Ref. [27]. This means that our 
model adequately reflects contributions of the extended 
interaction triads and thus should be useful in the fu- 
ture for a quantitative description of the transient region 
between turbulent spectra with the thermodynamic and 
the energy-flux equilibria. 



APPENDIX B: VELOCITY CIRCULATION AND 
POLARIZATION IN TURBULENCE 

In this Appendix we will calculate the velocity circu- 
lation r over a circular contour of radius R in classical 
turbulence with a power-law spectrum. Let the second 
order velocity correlation function (3D spectrum) in the 
k-space for isotropic homogeneous turbulence be 



Ft — C F 



vlk 



x-3 



k "■ 



(Bl) 



where C F = 27r 2 |a; — 3|, v T is the rms velocity in turbu- 
lence and fc* is a wavenumber of truncation from above 
for x < 3 (e.g. for the thermodynamic equilibrium with 
x = 0) and from below for x > 3 (e.g. for K41 turbulence 
with x = 11/3). Such a truncation is necessary for v T to 
be finite, and we will see below that the x = 3 boundary 
also separates different types of the scaling behavior of 
the velocity circulation. 
Consider the circulation 



/ uj n d 2 r, 
Jr 



(B2) 



where u> = V x v is vorticity, and the integral is taken 
around an arbitrary circle of radius R. Then 

(r 2 ) = J J K„ L0 2 , n ) d 2 n d 2 r 2 , (B3) 

Due to isotropy of the turbulence we may approximate 
(wi, n u; 2 , n ) = \{<jJi-u)2). With ri2 = ri - r 2 and 
{u) kl ■ u) k2 ) = 2k\ (27r) 3 <S(ki + k 2 ) F kl , we have 

<«i • o*> = //^r e^ k " 1+k ^ (u: kl ■ Wfca ) (B4) 
k . ri2 d 3 k _ C F v 2 T kr 3 /•° sin(fcn 2 ) 



2\k 2 F k e i 



k x ~3 



dk . 



(2wy 7r z n 2 J 

When 3 < x < 5, this integral converges and we have 

v 2 
ri2 



(wi • wa) = 2\x - 3\^(Kr 12 ) x - 3 \-g(4 - x) sin ^} , 

(B5) 



where Q{x) is the Gamma- function. Substituting this 
expression into (B3) and integrating we have 



(B6) 



where C x is an order-one constant (whose analytical de- 
pendence on index x is very complicated). 

For x < 3, the integral in (B4) diverges at the upper 
limit and, therefore, has to be truncated at the maximum 
wavenumber which, in this case, is fc*. We have, 



(u>i • u> 2 ) = 



C F vl {Kr 12 y 



3 rk,n 
JO 



■ sin y 



dy . (B7) 



The integral in this expression can be found in terms of 
the special functions whose asymptotical behavior can 
be readily obtained. This way one can show that the 
correlator (u>i ■ u) 2 ) decays in r\ 2 sufficiently fast, so that 
for k*R ^> 1 one can pass in the integral (B3) to the 
symmetric variables r + = |(ri+r 2 ) and r 12 = Y\— r 2 , use 
the polar coordinates and replace the upper integration 
limit for r\ 2 with infinity, 



<r 2 ) = 



2tt 2 R 2 



/ ( 

Jo 



oj\ ■ u) 2 ) r 12 dr\ 2 . 



Substituting (a>i • a> 2 ) from (B7) and integrating we have 



4tt 4 



v 2 R 2 . 



(B8) 



Interestingly, this expression is independent of both the 
spectrum exponent x and of the cutoff fc* . This fact has 
a simple physical interpretation. Suppose that the cor- 
relation length of vorticity field in turbulence d u is short 
so that l u <C R- Then the circle or radius R embraces 
N = R 2 /^ random and effectively independent vortex 
tubes each having radius ~ and circulation 7/ ~ u T f u . 
Because of the statistical independence of these tubes, 
(r 2 ) can be found using the Central Limit Theorem, sim- 
ilarly to the way we did it in the main text for the set of 
random quantized vortex lines. We have 

<r 2 )~ 7 2 A^ 2 i? 2 , 

which coincides, up to an order-one numerical factor, 
with expression (B8). Note that dependence on l u has 
dropped out, which corresponds to independence of ex- 
pression (B8) of x and fc*. 

For x > 5, the integral in (B4) diverges at the lower 
limit and, therefore, has to be truncated at the minimum 
wavenumber which, in this case, is again fc*. We have, 



(<*>i • w 2 ) 



7T 2 (X — 5) 



(B9) 



As we see, this correlator is independent of the distance 
r\ 2 which simply means that the correlation length in this 
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case is of order of the maximal lengthscale 1/fc* . The in- 
tegration in (B3) in this case reduces to the multiplication 
by the square of the circle area, i.e. 



vl hi 



■(vri? 2 ) 2 = 



2?r 2 (x- 3) 



vlk 2 R\ 



3 7T 2 (x-5) v ; 3(a:-5) 

(BIO) 

The R A scaling here coincides with the one obtained in 
the main text for a bundle of perfectly aligned (polarized) 
vortex lines. This is not surprising since the vorticity cor- 
relation length in this case ~ 1/ which is much greater 
than the contour size R. Another interesting effect to 
note is that (r 2 ) diverges for x — > 5. 

Let us now summarize all the cases. We have 



x < 3 (e.g. thermodynamic) : (T 2 ) 
3 < x < 5 (e.g. K41) : (r 2 ) 
x > 5 (e.g. smooth field) : (r 2 ) 



^, 2 i? 2 



2n 2 (x-3) 2 , 



x-1 



3(x- 5) 



hi R 4 



For polarization P (see definition in the main text) we 
have respectively 




x < 3 
3)/2, xe(3,5), 
x > 5 . 



(Bll) 



So, the thermal equilibrium state (x = 0) is not po- 
larized at all, P = 0, whilst Kolmogorov turbulence 
(x = 11/3) is partially polarized, P = 1/3. 



[1] R. J. Donnelly, Quantized Vortices in He II (Cambridge 

University Press, Cambridge, 1991) 
[2] Quantized Vortex Dynamics and Superfluid Turbulence, 

ed. by C. F. Barenghi et al., Lecture Notes in Physics 

571 (Springer- Verlag, Berlin, 2001) 
[3] W. F. Vinen and J. J. Niemela, J. Low Temp. Phys. 128, 

167 (2002) 

[4] S. N. Fisher, A. J. Hale, A. M. Guenault, and G. R. 

Pickett, Phys. Rev. Lett., 86, 244 (2001) 
[5] R. Hanninen, R. Blaauwgeers, V. B. Eltsov, A. P. Finnc, 

M. Krusius, E. V. Thuneberg, and G. E. Volovik, Phys. 

Rev. Lett. 90, 225301 (2003) 
[6] A. P. Finne, T. Araki, R. Blaauwgeers, V. B. Eltsov, N. 

B. Kopnin, M. Krusius, L. Skrbek, M. Tsubota, G. E. 

Volovik, Nature 424, 1022 (2003) 
[7] D. I. Bradley, D. O. Clubb, S. N. Fisher, A. M. Guenault, 

R. P. Haley, C. J. Matthews, G. R. Pickett, V. Tsepelin, 

and K. Zaki, Phys. Rev. Lett. 95, 035302 (2005) 
[8] V.B. Eltsov, A. P. Finne, R. Hanninen, J. Kopu, M. Kru- 
sius, M. Tsubota, and E. V. Thuneberg, Phys. Rev. Lett. 

96, 215302 (2006) 
[9] D.I. Bradley, D. O. Clubb, S. N. Fisher, A. M. Guenault, 

R. P. Haley, C. J. Matthews, G. R. Pickett, V. Tsepelin, 

and K. Zaki, Phys. Rev. Lett. 96, 035301 (2006) 
[10] T. Araki, M. Tsubota, and S. K. Nemirovskii, Phys. Rev. 

Letts. 89, 145301 (2002) 
[11] S.R. Stalp, J.J.Niemela, W.F.Vinen and R.J.Donnely, 

Physics of Fluids 14, Issue 4, pp. 1377-1379 (2002) 
[12] E.V. Kozik and B.V. Svistunov, Phys. Rev. Lett. 92, 

035301 (2004) 

[13] W. F. Vinen, M. Tsubota and A. Mitani, Phys. Rev. 



Lett. 91, 135301, (2003) 
[14] C. Connaughton and S. Nazarenko, Phys. Rev. Lett. 92, 

044501 (2004) 
[15] S. Nazarenko, JETP Lets., 84, 700 (2006) 
[16] B. V. Svistunov, Phys. Rev. B 52, 3647 (1995) 
[17] W. Thomson, Philos. Mag. 10, 155 (1880) 
[18] H. E. Hall, Proc. R. Soc. London, Ser. A 245, 546 (1958) 
[19] S.K. Nemirovskii, J. Low Temp. Phys. 142, 769 (2006). 
[20] K.W. Schwarz, Phys. Rev. B 31, 5782 (1985) and 38, 

2398 (1988) 

[21] H. Hasimoto, J. Fluid Mech, 51 pp 477-485(1972) 
[22] W. F. Vinen, Phys. Rev. B 61, 1410 (2000) 
[23] S. Nazarenko, JETP Letters. 84, issue 11, 700-702 
(2006). 

[24] S.A. Orzag, J. Fluid Mech. 41, 363 (1970) 
[25] W. J. T. Bos and J.-P. Bertoglio, Phys. of Fluids, 18, 
071701 (2006) 

[26] V.B. Eltsov, A. Golov, R. de Graaf, R. Hanninen, M. 
Krusius, V. L'vov, R.E. Solntsev, "Propagation of turbu- 
lent front in rotating superfluid" , International workshop 
on rotating superfluid turbulence, Propagation of turbu- 
lent front in rotating superfluid turbulence, Jerusalem, 
April, 2007. 

[27] V.I. Belinicher and V.S. L'vov, Zh. Eksp. Teor. Fiz., 93 
1269 (1987) 

[28] A. Eidelman, T. Elperin, A. Kapusta, N. Kleeorin, A. 
Krein, and I. Rogachevskii, Nonlinear Processes in Geo- 
physics 9 201 - 205 (2002) 

[29] E. V. Kozik and B. V. Svistunov, (2007), 
arXiv:cond-mat/0703047v2 [cond-mat. other] 



